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Abstract 

Context. The high number of planet discoveries made in the last years provides a good sample for statistical analysis, leading to some 
clues on the distributions of planet parameters, like masses and periods, at least in close proximity to the host star. We likely need to 
wait for the extremely large telescopes (ELTs) to have an overall view of the extrasolar planetary systems. Those facilities will finally 
ensure an overlap of the discovery space of direct and indirect techniques, which is desirable to completely understand the nature of 
the discovered objects, obtaining both orbital parameters and physical characterization. 

Aims. In this context it would be useful to have a tool that can be used for the interpretation of the present results obtained with various 
observing techniques, and also to predict what the outcomes would be of the future instruments. 

Methods. For this reason we built MESS: a Monte Carlo simulation code which uses either the results of the statistical analysis of the 
properties of discovered planets, or the results of the planet formation theories, to build synthetic planet populations fully described 
in terms of frequency, orbital elements and physical properties. They can then be used to either test the consistency of their properties 
with the observed population of planets given different detection techniques (radial velocity, imaging and astrometry) or to actually 
predict the expected number of planets for future surveys, as well as to optimize the future multi-techniques observations for their 
characterization down to telluric masses. 

Results. In addition to the code description, we present here some of its applications to actually probe the physical and orbital 
properties of a putative companion within the circumstellar disk of a given star and to test constrain the orbital distribution properties 
of a potential planet population around the members of the TW Hydrae association. Finally, using in its predictive mode, the synergy 
of future space and ground-based telescopes instrumentation has been investigated to identify the mass-period parameter space that 
will be probed in future surveys for giant and rocky planets 

Key words. Stars: brown dwarfs, planetary systems - Methods: data analysis, statistical 

1 . Introduction tential of the direct detections in the exploration of the outer re- 
gion of the planetary systems, also raising many question s about 

Many statistical studies have been done using information com- j^^^ ^^^j^ ^^-^^^^ ^^^j^ (^^^ lAbsil & Mawej |2009). This 

ing from more than a decade of extensive searches for exoplan- ^^^^^^ ^j^^ ^-^^^^ ^j^^ ^^^^ generation high co ntrast imaging in- 

ets, trying to answer questions either related to the properties of ^tr^^jnents like the Gemini Planet Imager (GPI: iMacintosh et alJ 

those obiects such as the mass, orbital period and eccentricity ,2007) and VLT/SPHERE (Spectro-Polarimetric High-contrast 

d Lineweaver&Gretheil |2003t | Cummmg et al. | |2008|), or about gxoplanet REsearch: Beuzit et al. 2008). These instruments will 

the relevance of the host star characteristics (mass, metallicity j^j^^j^ ^jj^^ ^^^^^^ ^^^^^ ^ systematic characterization to 

and binarity) on the final frequency and distabution of plan- ^^^j^^ 10 ^^jy j^^^ j^^j jii^itations (inner 

etary systems, (see Fischer & Valentj |2005|; I Santos etalj |2004t ^^^y.- ^„ ^^^^ ^^^^^^^ achievable), these instrument will 



| Johnsonetal. | |2007b . Since the most successful techniques (ra- ^^^^ gj^„t pj^^^t^ ^^^-^^ ^^^y ^^^^ ^^^^-^ ^^^^ 



dial velocity and transit) have focused on the inner (< 5 AU) ■ ^ for the ELTs facilities. A wide range of planetary 

environment of main sequence solai--type stars, most of the avail- ^^^^^^ separations, down to the rocky planets (and, in very 

able information on the frequency of planets concern this class f^yo^rable cases reaching the habitable zone), will be explored 

with 30-40 meter-class telescopes, finally allowing an overlap 

Recent discoveries of young distant planetary mass objects between the discovery spaces of direct and indirect techniques. 
with direct imaging (see e . g. Mar ois et alJ 12008: Kal as et alJ 
l2008t iLagrange et al.ll2008l) are giving us a first hint on the po- 

In this context it is useful and crucial to predict the perfor- 

* www.messthecode.com mances of the forthcoming instruments, not only in terms of 

** send offprint requests to bonavita@astro.utoronto.ca number of expected detections, but also trying to figure out what 
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will be the explored parameter space and even the possible syn- 
ergies between different discovery techniques. 

Here we present our Monte Carlo simulation code MESS, 
whose aim is to provide a flexible and reliable tool for the statis- 
tical analysis and prediction of the results of planet searches. 

It produces synthetic planet populations, deriving all the 
physical parameters of these planets together with the observ- 
ables that can be compared with the predicted capabilities of ex- 
isting or planned instruments. Such comparisons allow to derive 
subsets of fully characterized detectable planets, as well as a 
snapshot of what the evolution of the sample of detected planets 
would be in the next years. 

A detailed description of the code, and of all the assumptions 
which constitute its basis, is given in Sect. |2] while in Sect. [3] 
we present the different operation modes of the code and their 
applications. Although the MESS has been built, and it has been 
so far applied, only to analyze and/or predict the results of direct 
imaging surveys, an extension of the code to different techniques 
is planned. The first attempt in this direction are presented in 
Sec. 13.31 Conclusions and suggestions for further work will be 
finally drawn in Sect. |4] 



MESS (Multi-purpose Exoplanet Simulation 
System) 



Over the past years , several groti ps ("Kasper et alj 
Lafreniere etaP 120081: IChauvin et aP 12010: Nielsen & Close! 



2007t 



2O10ir initiated statistical analysis to constrain the physical and 



orbital properties (mass, period, eccentricity distributions) of the 
giant planet population at large separations. They developed sta- 
tistical analysis tools appropriate to exploit the performances of 
deep imaging surveys. They also tested the consistency of var- 
ious sets of parametric distributions of planet properties, using 
the specific case of a null detection. The first assumption of 
these tools is that planet mass and period distributions coming 
from the statistical results of RV studies at short per iod (see e.g. 
iLineweaver & Gretheii 120031 : ICumming et al.ll2008l) can be ex- 
trapolated and normalized to obtain information on more distant 
planets. Despite the model-dependency on the mass predictions, 
the approach is attractive for exploiting the complete set of de- 
tection performances of the survey and characterizing the outer 
portions of exo-planetary systems. 

With all of this in mind, we tried to go a step further, cre- 
ating a Multi-purpose Exo-planet Simulation System (hereafter 
MESS) to be applied also to other techniques than direct imag- 
ing, also using the information coming from the planetary for- 
mation theories. 

The code is written in IDL and can be downloaded from 
www.messthecode.com 

The basics operations performed by the code are the follow- 
ing: 

1 . it generates a synthetic population of planets, including all 
the orbital elements, either using the planet mass and period 
distributions coming from the statistical results of RV studies 
or the outcome of the planetary formation theories. 

2. taking into account the characteristics of the host star and 
of the planetary orbit, it calculates all the observable quan- 
tities needed for the comparison with the instrument perfor- 
mances, such as radial velocity (RV) and astrometric signal, 
planet/star contrast, degree of polarization, etc. 

3. given the detection capability relation of an instrument, ei- 
ther already available or planned, it selects a sub-sample of 



fully characterized detectable planets, which characteristics 
can then be analyzed. 

The code then assumes a given star population, a planet popu- 
lation with associated physical and orbital properties based on 
a theoretical or semi-empirical approach, the corresponding ob- 
servables for different observing techniques, finally generate a 
synthetic population of planets to be compared with the instru- 
mental detection performances. Each step is described hereafter. 

2.1. Star population 

The first input of the MESS is a sample of A^star stars, which have 
been targeted for planet searches or which are part of a sample 
for future observations. Various stellar parameters are assumed 
to be known, such as the apparent magnitude, the distance, the 
luminosity, the spectral type, the mass, the age, the metallicity, 
etc. 

Fig. [1] shows the characteristics of a sample of 600 nearby 
(d < 20 pc) stars se lected from the Hipparcos catalogue 
dPerrvman & ESAlll997h and used to build the synthetic pop- 
ulation showed in Fig. [3] 
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Figure 1. Principal ciiaracteristics of the sample of nearby stars used 
to built the example synthetic population. Upper Left: Apparent 
Magnitude in the J Band vs distance in pc. Upper Right: Stellar Mass 
(Mo) vs distance in pc. Lower Left: Histogram of stellar ages (Myrs). 
Lower Right: Histogram of stellar masses (Mq). 



2.1 .1 . Binarity module 

MESS also gives the possibility of taking into account the pres- 
ence of one (or more) additional stellar companions, in the anal- 
ysis. If a star in the sample is flagged as binary, the code uses the 
information about the binary orbit (if available) to compute the 
critical semi-major axis for the dynamical stability of the system. 
This set the limiting value that the semi-major axis of a planet 
can attain and still maintain its orbital stability, as a function of 
the mass-ratio and orbital elements of the binary, as shown by 
iHolman & Wiegertl(ll999l) . 

Both the case of circumstellar (or Satellite S-type) and cir- 
cumbinary (or Planet P-type) orbit are considered, and the criti- 
cal s emi-major axis is compute d using Eq.[T]and|2]respectively, 
from lHolman & Wiegerll (Il999l) . 
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adai, = 0.464- 0.38 /i- 0.361 ei-H 0.586 yuei 

+ 0. 150 yu^ -0.198 jue^ (1) 

ajat - 1.60 -H 4.12/i + 5.10 gfe - 4.27 yu e/, -H 

- 5.09 - 2.22 el + 4.61e^ (2) 

In both the equations, Uc is critical semi-major axis , yU - 
Ml /(Ml + M2), Ob and ei, are the semi-major axis and eccen- 
tricity of the binary, and Mi and M2 are the masses of the pri- 
mary and secondary stars, respectively. If not available from 
literature, the eccentricity is assumed to be eh - 0.36, re- 
ported a s mean value for the eccen tricity of a binary sys- 
tem, by iDuquennov & Mavoil (Il99lh . In case the value of 
the semi-major axis is not available, then the code estimates 
it as Ub = 1.31 p (arcse c) d (pc) (see iFischer et alj l2002t 
iDuauennov & Mavoilll99llPl 

Note that in the first case (S-type orbit), set the maxi- 
mum value that the semi-major axis of a planet can assume, 
before compromising the stability, while it represents the min- 
imum value of the semi-major axis of a stable planet, in the case 
of a P-type orbit. 



2.2. Planet population 

The core of the code is the generation of the synthetic planets, 
that are fully characterized, both in terms of orbital parame- 
ter, and physical characteristics. Depending on the goal of the 
study one can choose between a Semi-empirical approach or 
a Theoretical Approach. These different approaches makes the 
code suitable to constraint the planet properties under different 
assumptions, but also to test model predictions. 
If the Theoretical Approach is chosen, masses and period val- 
ues selected from a synthetic population prov ided by the out- 
put of the planetary formation models (see e.g. M ordasini et alj 
[2009) are given as input. In this case all the orbital characteris- 
tics are also provided, together with the physical properties of 
each planet, so no random generation is needed, and the code 
only evaluates the observable and compares them to the provided 
detectability relations. Different populations of planets obtained 
assuming different stellar masses and metallicity values can be 
selected according with the characteristics of the real star in the 
sample, to take into account the effects of the stellar characteris- 
tics on the planet formation processes^. 

The Semi-empirical Approach uses the power law distribu- 
tions in Eq. [3] and |4] for the mass and semi-major axis of the 
planets as retrieved from the statistical analysis of the properties 
of the planets discovered so far to generate a seed population of 
Nseed values of masses and periods (see Sec. 12.2.11 . 



dN 
d{Mp) 
dN_ 

Ip 



(3) 



(4) 



' Note that lDuguennov & Mavoj ( Il99lh refers to solar-type star mul- 
tiplicity 

^ The results discussed in this paper has been obtained using mainly 
the Semi-empirical Approach. An extensive use of the Theoretical 
Approach, using as input the newest outcomes of the Bern formation 
models (Mordasini et al. 2011), will be the subject of a forthcoming 
paper. 



The user can also set a pre-determined grid of masses- 
periods and feed it to the code, without any assumption on the 
distributions. This would be the case if, for example, bound- 
aries on the mass/semi-major axis space where planets can form 
are to be set using the outcomes of a formation model(see e.g 
iMordasini et al.ll2010h . excluding from the sample the planets 
not compatible with the theory. 

If the Semi-empirical approach is used, mass, orbital param- 
eters, as well as temperature and radius of the planets, are ob- 
tained based on the assumption described in the next sections 

2.2.1. Mass-period seed generation 

If the semi-empirical approach is chosen, the power law distribu- 
tions are fed to the Monte-Carlo core of the code, that randomly 
generates a fixed number of mass-period pairs. Both the plan- 
etary mass and period ranges can be given as inputs, together 
with the power-law exponents. In a typical setup, the power-law 
exponents a re assumed to be a = - 1.31 and j8 - -0.74 respec- 
tively, from lCumming et aP (12008 '). The planetary masses span 
the range between 0.6 M Earth and 15 My,,,,, and the periods (P) 
are chosen between 2.5 days and 350 years (corresponding to 50 
AU for 1 Mq star). 

A scaling of the planetary mass, and even of the period, with 
the stellar mass can be als o introduced, according to recent re- 
sults (e.g. iLovis & Mavoii l2b07). In addition, a dependence of 
the planet f requency on the s tellar metallicity may also be con- 
sidered (see lFischer & Valentiii2005i) . 

2.2.2. Evaluation of tine orbital parameters 

For each mass-period pair in the seed generation, the code eval- 
uates the semi-major axis computed using Kepler's third law, 
using the mass of each star in the input sample. Then, it gen- 
erates Ngen values of all the orbital parameters: eccentricity (e), 
inclination (/), longitude of periastron (w), longitude of ascend- 
ing node (Q), and time of periastron passage (To). By default, all 
these parameters are randomly generated following an uniform 
distribution 0. The eccentricity distribution is cut al e - 0.6 as 
suggested by the results of the RV surveys (see Gumming et al. 
2008). This also allows to control possible bias towards high ec- 
centricity planets that could affect the results of Direct Imaging 
surveys. A full discussion of the impact of the eccentricity dis- 
tribution on the simulations results is held in Sec. 3.4. 

The date of observation is also required. If not available from 
the real data, an epoch of observation, tobs, is generated over a 
time-span chosen according with the considered instrument. 

The code also offers the possibility to fix each orbital pa- 
rameter to known or predicted values, for all the planets in the 
population. 

The coordinates, x and y, of the projected orbit on the plane 
perpendicular to the li ne of sight, are finally computed using the 
ephemeris formulae of lHeinta(ll978l) . reported in Eq.|5]to|7] 



X ^ AX + FY 
y = BX-\-GY 
X - cos E - e 
Y = Vl -e2 sinE 



P 



^Jx^Ty^ 



(5) 
(6) 

(7) 



^ Note that in the case of the inclination, cos / and not i itself is uni- 
formly generated by the code 
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where X and Y are the coordinates of the orbit (Eq. |6ll, p is 
the projected separation, and A, B, F, G are the Thiele-Innes ele- 
ments, which can be obtained from the classical ones (the semi- 
major axis a, o), Q., and /) using Eq. |8] 

A - a (cos (jj cos Q - sin <jj sin Q cos /) 

B - a (cos d) sin Q + sin cj cos Q cos /) (8) 

F - a(-sm(L) cos Q - cos co sin O cos 

G - a (- sin sin Q H- cos w cos O cos • 

In these equations, E is the eccentric anomaly (obtained from the 
mean anomaly M ( Eq. |9ll using Eq. [TOl l and v the true anomaly 
(Eq.in): 

M = | ^°'^~^° j27r (9) 

El) - M + esmM + — sm 2M 
2 

Mq = Eq - e sin £0 
£ = £o + (M-Mo))/(l -ecos£o) (10) 
tanv/2 = V(l +e)/(l -e)tan£/2 (11) 

The projected separation, p (in arcsec), can be obtained ei- 
ther using Eq. Qor Eq. [12] (which gives also an estimate of the 
radius vector: r), then dividing for the star distance. 

p = r cos (v -H w) sec (0 - Q) (12) 
r = fl(l -e^)/(l +ecosv) 




1 10 100 1000 
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Figure 2. Summary of the planetary Mass - Radius relations adopted 
for the different mass ranges. All the model computation are made as- 
suming a host star of 1 Mq, and the semi-major axis value is fixed to 
5 AU. Filled symbols corresponds to known transiting planets; open 
symbols are for Solar System planets. 

2. Equations [TSjandfTSlfrom'Fortne v et alj (l2007h are used for 
the smallest planets (M < lOMEaith)- These are either: 

R = (0.0912 /m/ + 0.1603)(logM)2 
-i-(0.3330 /m/-H 0.7387) log M 
H-(0.4639 /m/+ 1.1193) (15) 



2.2.3. Planet Temperature 

Since we aim at consider both the thermal and reflected flux of 
the planets, we need two different estimates of the temperature. 
The first one is the internal temperature. Tint, coming from the 
evolutionary models (see e.g. Baraffe et al.. 2003i) . The second 
one is the equilibri um temperature, Teq, obtained through Eq.[T3] 
(from lSudarskv et al. 2003.) 



- eq 



(1-Ab)L, 



I6ncra^ 



(13) 



where L« is the star luminosity. The Bond alb edo Ab is assumed 
to be 0.35 in the J band (Jupiter value, see iHanel et al.|[T981h 
and it is randomly generated between 0.3 and 0. 52 in the visible 
(the l atter being the Jupiter albedo in V band.see lSudarskv et alj 
I200I 

Our final assumed value for the effective temperature of the 
planet T^ff is given by: 



7^4 _ 7^4 
^ eff ~ ^ int 



+ T1 



(14) 



2.2.4. Planet radius 



MESS uses the approach developed bv lFortnev et alJ (l2007h to 
evaluate the planetary radius. Practically the radius is assumed 
to depend on the planet mass, with the following recipes: 

1. For Jupiter-like planets (M > lOOMearth ), an interpo- 
lation is performed within the published values given by 
iFortnev et al.l (l2007h . Values of age and distance of each 
star are entered, yielding a value for Rcas- A core mass of 
lOMEarth is assumed. 



R = (0.0592 rOT/ + 0.0975)(logM)2 
-l-(0.2337 nnf + 0.4938) log M 
H-(0.3102 nnf + 0.7932) 



(16) 



for ice/rock and rock/iron planets, respectively. In these 
equations, R is in /?Earth and M is in MEarth, while imf is 
the ice mass fraction (1.0 for pure ice and 0.0 for pure rock) 
and nnf is the rock mass fraction (1.0 for pure rock and 0.0 
for pure iron). In the typical MESS setup, the ice/rocky or 
rocky/iron fraction is set to 0.3 (50% of chance for a planet 
being mainly icy or rocky). 
3. Finally, predictions are uncertain for the Neptune-like plan- 
ets, where the transition between the two relations described 
above should occur. The most sensible approach seems to be 
to fit the mass-radius relation of the Solar System in the same 
mass-range (10-40MEai-th )■ This procedure provides a good 
agreement with the radii of Uranus and Neptune and of the 
few transiting Neptunes confirmed so far (as listed by The 
Extrasolar Planet EncyclopaedicQ see Fig. |2]l. 

The resulting mass-radius relations are showed in Fig. |2] 
with over-plotted the data corresponding to the planets discov- 
ered with the transit technique and the planets from our Solar 
System, for comparison. 



2.3. Predicted observables 

Having in hands the full set of orbital and physical parameters 
of the planets, the code then provides an estimate of observable 



www.exoplanet.eu 



Bonavita et al.: MESS (Multi-purpose Exoplanet Simulation System) 



5 



quantities such as the luminosity contrast or the degree of po- 
larization, needed for direct observations, but also quantifies the 
indirect effects of the presence of the planet, providing a measure 
of the semi-amplitude of radial velocity (RV) and the astrometric 
signal. 

2.3.1. Planet/Star contrast 

MESS gives an estimate of both the intrinsic and reflected flux, 
in the selected band, for each planet. Throughout the paper we 
will refer to the planets which luminosity is dominated by the 
intrinsic contribution as self-luminous or warm planets, as op- 
posed to the cold planets for which the reflected light provides 
most of the contribution to the planet/star contrast. 

The intrinsic emission is estimated using the prediction of 
evolutionary models at the age of the star (assumed to be also 
the age of the system). To this purpose two classes of models 
can be considered, based on different assumptions on th e initial 
conditions: Hot Start mod els (Chabrieret al. 2000; Baraf fe et alj 
l2003t [S aumon & M arlev 2008), which consider an initial spher- 
ical contracting state; and Core Accretion models ( Marl ev et alJ 
120071; iFortnev et al.ll2008l) . which couple planetary thermal evo- 
lution to the predicted core mass and thermal structure of a core- 
accretion planet formation model. 

In the following, we only consider the results obtained using 
the hot start models for the nearby sample. However, the problem 
of the initial condition and the uncertainties on the stellar ages 
are among the main limitations, in case of young stellar sam- 
ples, not only for our code, but also for a ny kind of study that 
uses the same kind of approach (see e.g. IChauvin et al.ll2010t 
iBonavita et al.ii2010 ,for a detailed discussion). These limitations 
also apply to the theoretical approach, if the evolutionary mod- 
els are used to evaluate the planet intrinsic luminosity and radii 
produced by the models, as in Mordasini et al. (2010, 2011). 

For the evaluation of the reflected light, we scaled the Jupiter 
value, according with the planet radius (expressed in Jupiter 
radii), semi-major axis, albedo and illuminated fraction of the 
planet. This last contribution is computed through a phase de- 
pendent term, 0(^6), which is given by Eq.[T7](see Brown 2004), 
where /3 is the phase angle (angle at companion between star and 
the observer) and z = r sin (v + w) is the radial coordinate of the 
radius vector. 

(b(/3) = [sm/3 + (7T+/3)cos/3]/7T (17) 

The Jupiter/Sun contrast is obtained using Eq. [18] which 
gives an estimate of the fraction of stellar light captured by a 
planet, depending on the values of the planet radius, semi-major 
axis and geometrical albedo, being <t(/?) - 1 (at opposition). 

(Ljup/i.)Ref = Aj„p^ = 2.5 X 10-' (18) 

Jup 

Where A jup - 0.35 is the valu e of the Jupiter albedo in the the 
J-Band, (see lHanel et aPflOSl . 

Then we end with a final value of the contrast in reflected 
light given by EafT9l 

(Lp/L.).,.(L.p/L.)^^^0(^)^g;^ (19) 

As a consequence of Eq. 19, the results of MESS will be sensi- 
tive to the choice of A^, especially for the cold planets, in which 
the contribution of the reflected light is dominant. Following the 



out comes of Jupite r observations and theoretical models (See 
e.g. lBurrowsll2004h . we decided to uniformly generate the val- 
ues of the albedo between 0.2 and 0.7. The code anyways offers 
the option to fix the value of the albedo to a chosen value, for all 
the planets in the generation. A test of the impact of the choice 
of the albedo value on the redults of the simulations is presented 
in Sec. 3.4. 

2.3.2. RV and astrometric signal 

The indirect effects of the presence of the planet, such as the 
semi-amplitude of radial velocity (RV) variations and the astro- 
metric signal can be inferred, knowing all the orbital character- 
istics for each planet. 

2.3.3. Degree of polarization 

The degree of polarization 11 is assumed to be of the form (see 
e.g.fStenflo 2005): 

n = nmaxX(l-cos^/3)/(l+cos2/3) (20) 

where Dmax is the maximum polarization value (which is as- 
sumed to be randomly generated between 0.1 and 0.3), and fi 
is the same as in Eq. [17] Then the contrast due to the polarized 
light of the planets is 11 times the contribution in reflected light 
evaluated with Eq. [19] 

2.4. Planet population synthesis 

Depending on the purpose of the analysis, the code can generate 
the planet population in two different ways: 

a) Full population: the value of Nseed sets the spacing of the 
mass-period grid, and for each point on it Ngen planets are 
generated, ending with NseedXNge,, planets per each star. The 
population for each star is saved in an independent file. This 
approach is useful for the statistical analysis of existing data, 
since in this case MESS provides the fraction of detectable 
planets per star, which can be used to derive the global prob- 
ability of finding a planet over the whole target list. This can 
be then compared with the real results. 

b) Reduced population: only one orbit is generated for each 
point in the mass-period grid. Nge„ in this case sets the num- 
ber of planet in a planetary system associated with each stail^. 
The final population is then composed by Nstar x Nseed plan- 
ets, and all the planets are saved together in one file. Then the 
predicted detection performances of a given instrument can 
be used, to derive the population of objects that are expected 
to be detected around each star, if the whole input sample is 
observed. 

As an example, we generated a reduced population (assum- 
ing 5 planets per star) of planets around the stars of the nearby 
sample described in Sec. 12. II We choose the semi-emp irical _a p- 
proach, and used the typical setup we discussed in Sec. 12.2. lln 

Fig Sshows the position of the planets in the mass vs semi- 
major axis plane. 

^ note that no consideration on the planet stability is made, and to the 
purpose of the analysis each planet is considered separately 

* Note that the whole calculation of the physical characteristics and 
observables described in Sec. 2.2.3 to 2.3.2 can be skipped (with con- 
siderable gain in computing speed), the code providing in this case only 
the orbital elements 
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Figure 3. Mass semi-major axis distribution of the synthetic planets in 
the populations generated by MESS using the semi-empirical approach. 
The different classes of planets (see text) are plotted using different 
colours: red/orange for the warm/cold Jupiters, green for the Neptune 
like planets, blue for the rocky planets. 
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Figure 4. Distribution of radial velocity vs. period of the synthetic 
planets for the population showed in Fig. [3] 

The planets are separated into the three classes, using differ- 
ent colours: 

- Giant (or Jupiter-like) planets (Mpianet > 40MEarth)- A 
distinction between Cold Jupiters (orange dots) and Warm 
Jupiters (red dots), as defined in Sec. 12.3.11 is also made. 

- Neptune-like planets (lOMEaith < A^pianet < 40MEaith: green 
dots) 

- Rocky planets (Mpianet < lOMEarth: blue dots) 

The distribution of the observable quantities for the planet 
showed in Fig. [3] are summarized in Fig.|4]|5]and|6] 

2.5. Instrument detection performances 

The last step is the comparison of the observables of the gen- 
erated synthetic planets with the detection limits of different 
observing techniques, with the possibility to actually combine 
them. It is important, especially in case of comparative studies. 



Figure 5. Distribution of the astrometric signal vs. period of the syn- 
thetic planets for the population showed in Fig.|3] 

10"* I ' ' — ' — ' ' — ' — ' — I 

Worm Jupiters • 
• Cold Jupiters • 

» * • Neptunes • 




projected seporotlon (orcsec) 



Figure 6. Planet/star contrast vs. projected separation of the same plan- 
ets showed in Fig. [3] 

and since the MESS does not produce the detection limits, to 
make sure that the detection performances that are fed as input 
to the code have been estimated by correctly taking into account 
each instrumental biases, specific to each technique, and the stel- 
lar characteristics. In the context of the MESS applications, the 
code has been extensively used considering two possible inputs 
for the detection performances: 

- The ID mode, which selects the detectable planets using a 
threshold or a curve giving the lower detection limits (RV, as- 
trometric precisions or contrast performances) as a function 
of the period, the semi-major axis, the angular separations 
etc., defined by the instrumental capabilities 

- The 2D mode, which is especially built for the analysis of the 
performances of the Deep Imaging instruments. This mode 
takes advantage from the knowledge of all the orbital ele- 
ments of the planets, to place them on a two dimensional 
detection map. This mode allows using all the spatial infor- 
mation stored in the images. Using the whole 2D map not 
only allows to take into account possible peculiar character- 
istics of the circumstellar environment, such as the presence 
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of disks, but also prevent to under/overestimate the contrast 
curve depending on the method chosen for the extraction it- 
self (see feonavita et aL,2010i) . 

3. Applications 

Once the synthetic population of planets has been created, the 
next step is to compare the characteristics of the generated plan- 
ets with the detection limits appropriate for the instrument under 
consideration. 

MESS offers three different operation modes (OM), depend- 
ing on which kind of analysis is needed. 

1. The Hybrid Mode (MESS.HM) which is the most flexible 
one, and an be used to probe the physical and orbital proper- 
ties of a putative companion around one given system based 
on the combination of different techniques, a priori informa- 
tion on the possible orbit given the presence of other planets 
or circumstellar disk. 

2. The Statistical Analysis Mode (MESS.SAM), which is 
built for the analysis of real data and uses the full popula- 
tion defined in Sect. 12.41 It enables to test different set of 
planet populations or constrain the maximum occurrence of 
planets for a given population that would be consistent with 
the results of detection and/or null-detection of a complete 
survey of a large target sample. 

3. The Predictive Mode (MESS_PM), which starts from the 
reduced population (see Sect. 12.4b . and given the predicted 
performances of a planned instruments, can be used to select 
the most suitable targets given the science goals of the instru- 
ment itself, to test the results of different observing strategies 
and finally to foresee possible synergies with other instru- 
ments. 

3.1. Single object cliaracterization 

The first and more versatile MESS mode is the so-called Hybrid 
mode. This mode can be used for the study of particularly inter- 
esting targets, or to test specific hypothesis. It allows for example 
to take into account all available informations about the orbit of a 
planet already discovered around the target, in order to put con- 
straints on the planet generation. A preliminary version of this 
mode has been used to put constraint on the presence of a plan- 
etary companion embedded in the disk surrounding the T-Tauri 
star LkCal5 (see Bonavita et al. 2010). 

We present here an analogous analysis made for TWA 1 1 . 
This star has been fou nd to be surr ounded by a debris dis k by 
[Schneider et an(l2009h . Using STIS. ISchneider et al.l(l2009l) pro- 
vided a full characterization of the disk geometry, and suggested 
a possible unseen companion responsible for some of the ob- 
served properties. We then decided to use MESS_HM to verify 
which kind of constraints can be put using the VLT-NACO ob- 
servations of this star 

A pixel-to-pixel 2D noise map was estimated from the re- 
duced NACO images, using a sliding box of 5 x 5 pixels over the 
whole FoV. We then considered a 6 cr threshold to build the final 
detection limit maps to be used for the statistical analysis. These 
maps were also converted in terms of minimum mass map us- 
ing the evolutionary model predictions at the age of the_system 
Fig. 7 shows an example of the resulting sensitivity maj 
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Figure 7. 2D map giving the values of the minimum mass of de- 
tectable companions (6cr) as a function of projected separation, 
around TWA 1 1 



We considered only circular orbits coplanar with the disk, 
with an inclination and a longitude of the ascending node fixed 
by the disk properties reported in Schneider et al., (200^): ipisk = 
75.88 + 0.16, Q = PA±9Q = (21. 1 + 90). ISchneider et al.1 (l2009h 
also estimated the inner and outer boundary of the disk to be 
0.515 and 2.114 arc sees respectively, cor responding to 37 an d 
154 AU at the target distance (72.8 pc, see I Van Leeuwenl2007h . 
TWA 11 is also k nown to have a ste llar companion at p = 7.7" 
(lJuraet al.ll 19951) . As pointed out bv lSchneider et al.1 (l2009h . the 
value of the outer boundary of the disk is consistent with the 
presence of the companion. Using Eq. [T] we in fact obtained 
a value for the critical semi-major axis for the planet stability 
{acrit) of about 165 AU. Taking into account these constraints, 
we set the range of explored semi-major axes to 35-160 AU. 
The results of our simulations, in terms of non detection proba- 
bility maps as a function of the companion mass and semi-major 
axes, are sh own in Fig. 8. The disk boundaries are also shown, as 
reported bv lSchneider et al.1 (l2009h . 

Is it clear that with the NACO images we are not able to put 
strong constraint on planetary-mass objects, but surely low-star 
companions and brown dwarfs more massive than 3QMjup can 
be excluded at a > 35 AU and IQMju,, ones for a > \QQAU. 



' Note that the decreasing values of the non-detection probability at 
separations lower than 30-50 AU are due to systematic errors. In fact 
the detection limit drops to unrealistic low values really close to the star 



3.2. Statistical Analysis of a Survey 

3.2.1. Testing the planet population assumption 

The MESS-SAM operation mode allows to test the consistency 
of various sets of (mass, eccentricity, semi-major axes) paramet- 
ric distributions of a planet population with observational data. 
Given the detection performances of a survey, the frequency of 
detected simulated planets (over the complete sample) enables 
derivation of the probability of non-detection of a given planet 
population associated with a normalized distribution set. Then 
the comparison with the survey results tests directly the disagree- 
ment with observations at an appropriate level of confidence. 

As an example of the use of SAM@MES, statistical analysis 
mode, we present the analysis of a small sample of young neigh- 
bourhood stars, part of the TWA association, and observed with 
NACO@VLT. These stars are part of a bigger sample for which 
the observations and statistical analysis, done with a preliminary 
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Figure 8. Non detection probability map of a faint companion 
around TWAll as a function of its mass and semi-major axis, 
in the case of a circular orbit. Inclination and longitude of the 
ascending node have been fixed using the disk properties ioi^k - 
75.88 + 0.16, = PA + 90 = (27.1 + 90). 



version of MESS, have been presented bv lChauvin et al. 
The characteristics of the stars in the sample are listed in Tab. 1 . 
We present a new analysis of these target, done with the 2D mod- 
ule ofMESS.SAM. 

2D minimum mass maps were obtained with the same 
method used for the analysis of TWA 1 1 (see Sec. 13.1b for all 
the stars in the sample. We then used MESS_SAM to calculate 
the detection probability (Pd) of companions of various masses 
and orbital parameters (semi-major axis a, eccentricities e, incli- 
nation /, longitude of the ascending node Q., longitude of perias- 
tron o) and time of periastron passage Tp). We used the empirical 
approach, generating a full population of 10.000 planets for each 
target, with a mass range spanning between 0.3 and 30 Mj,ip and 
a cut-off in semi-major axis of 100 AU. 

Each simulated companion was placed on the 2D minimum 
mass map according to its position on the projected orbit to test 
its detectability, comparing its mass with the minimum value 
achievable at the same position in the FoV. 

Only circumbinary planets were considered around TWA 22, 
adopting the total mass of the system as Mstar-Jn fac t the binary 
being so close (p - 0.1" see lBonnefovetal]|2009h leads to a 
value of the critical semi-major axis for circumstellar planets acs 
of only 0.456 AU and of 8.395 AU for the circumbinary ones. 

Two sets of indices for the power-law distribution were 
tested: 



1. The ones derived bv ICummin^ et alJ (l2008h (CM08): a = 
-1.31, j6 = -0.74 

2. The ones derived bv lLineweaver & Grethed (l2003l) (LW03): 
a = -1.81, jS = -0.30 

Finally, fixing a and to the CM08 values, we also intro- 
duced different values for the scaling of the planetary mass with 
the primary mass. 

The results of all these simulations are summarized in Fig. 9. 

3.2.2. Estimate of the frequency of giant planets 

A second more general use is to constrain the exoplanet fraction 
/ within the physical separation and mass probed by a survey. 
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Figure 9. Non-detection probability for the stars listed in Tab. 1 , 
based on various sets of period and mass distribution. Mass and 
period distribution are extrapolated and normalized from RV 
studies. Top: Variation of the non detection probability using 
two different sets of power-law distributions (see text). Bottom: 
variation o f the non-detection pro bability fixing a - -1.31 and 
j3 = -0.74 ('Cumming.etani2008) and different scaling the mass 
of the planet with the primary mass. 



in the case of null or positive detections. Contrary to what was 
assumed before, / becomes an output of the simulation, which 
actually depends on the assumed (mass, period, eccentricity) dis- 
tributions of the giant planet population. This statistical analysis 
aims at determining /, within a confidence range, as a function 
of mass and semi-major axis, given a set of individual detection 
probabilities pj directly linked to the detection limits of each star 
observed during the survey and to the considered giant planet 
distributions. 

The probability of planet detection for a survey of stars 
can in fact be described by a binomial distribution, given a suc- 
cess probability fpj, with / being the fraction of stars with plan- 
ets, and Pj the individual detection probabilities of detecting a 
planet if present around the star j. Each individual pj can be 
replaced by {pj}, the mean survey detection probability of de- 
tecting a planet if present. Finally, assuming that the number of 
expected detected planets is small compared to the number of 
stars observed (fipj) « 1), the binomial distribution can be 
approximated by a Poisson distribution to derive a simple ana- 
lytical solution for the exoplanet fraction upper limit /max for a 
given level of confidence CL: 



•ln(l - CL) 

N{pj) 



(21) 



Fig. 10 shows the results obtained applying this module at 
the sample of stars listed in Tab. 1 . 

Although the significance of our results is not really high, 
given the small size of the sample, they are stil l in agreement 
with th e results of the whole analysis presented bv lChauvin et al.l 
(1201 oh and with the results of the other deep imaging surveys 
(see e.g.^ielsen & Close 2010; Lafreniere et al. 2007)0 



* An extensive analysis, with MESS.SAM, of the results of the major 
deep imaging surveys published in the last decade is ongoing, and will 
be presented in a forthcoming dedicated paper. 
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Table 1. Sample of TWA stars considered in our analysis. In addition to name, coordinates, galactic latitude (b), spectral type, 
distance, V and K photometry, the observing mode direct imaging (DI) or coronagraphy (COR) and the status of the primary 
(single, binary Bin, triple) are also listed 



Name 


a 


6 


SpT 


Mass 


d 


Age 


V 


K 


Mode 


Notes 






[J2000] 


[J2000] 




Mo 


(pc) 


(Myr) 


(mag) 


(mag) 








TWA22 


10 17 26.9 


-53 54 28 


M5 


0.15 


18 


8 


13.2 


7.69 


DI/CI, S27, Ks 


Bin (p = 


0.1") 


TWA 14 


11 13 26.3 


-45 23 43 


MO 


0.55 


63 


8 


13.8 


8.50 


DI/CI, S27, Ks 






TWA12 


11 21 05.6 


-38 45 16 


M2 


0.30 


32 


8 


13.6 


8.05 


DI/CI, S27, Ks 






TWA 19 


1 1 47 24.6 


-49 53 03 


G5 


1.50 


104 


8 


9.1 


7.51 


DI/CI, SI 3, H 






TWA23 


12 07 27.4 


-32 47 


Ml 


0.40 


37 


8 


12.7 


7.75 


DI/CI, SI 3, H 






Twa25 


12 15 30.7 


-39 48 42 


M5 


0.15 


44 


8 


11.4 


7.31 


DI/CI, S27, Ks 






TWA 11 


12 36 01.0 


-39 52 10 


AO 


2.10 


67 


8 


5.8 


5.77 


DI/CI, S27, H 


Bin(p = 


7.7"), Star with disk 


TwalV 


13 20 45.4 


-46 11 38 


K5 


1.00 
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9.01 


DI/CI, SI 3, H 
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Figure 10. Top: Survey mean detection probability derived as 
a function of the semi-major axis , assuming para metric mass 
and period distributions derived bv lCumming et al.l (i2008) . The 
results are reported for individual masses: 1.5,3,6,9,12,15,30 
Mjup. The integrated probability for the planetary mass regime 
is shown with the thick green line. Bottom: Planet fraction upper 
limit derived as a function of semi-major axis, given the same 
mass and period distributions. 




Semi-major axis (AU) 

Figure 11. Non detection probability map, for 41 Ari (HIP 
13209). 



sive d isks. The full analy sis is presented in detail in the survey 
paper dJanson et al.l201 ih . An example of a detection probability 
map in mass versus semi-major axis space is shown in Fig. 11. 



3.2.3. Theoretical approach 

The MESS_SAM can also be used to test the predictions of spe- 
cific planet formation theories. An extensive use of this OM has 
been made to analyze a sample of massive stars (B-type and 
early A-type) observed with NIRI, to test the applicability of 
planet formation by disk instability in those systems. Starting 
from a uniform mass versus semi-major axis grid with a sam- 
pling of 5 AU in semi-major axis and 1 Mjup in mass, lO'* or- 
bits w ere generated for each grid point. M odels of disc insta- 
biUty dBell et alJll997l:lMordasini et alj|20"l lj,and Klahr et al., in 
prep.) were then used to provide boundaries in the mass ver- 
sus semi-major axis space, within which sub-stellar companions 
can form by this mechanism. These boundaries were dependent 
on the stellar properties, and so appropriate values should be 
used for each target in the sample. The planets falling within 
the allowed range were subsequently evaluated against the ID 
detection limits from the high-contrast images of the survey. In 
this way, by testing a range of planet distributions within the set 
boundaries, meaningful limits could be placed on the frequency 
of planet and brown dwarf formation by disk instability in mas- 



3.3. Predictive mode 

Beside the analysis of the real data, MESS can be also used to 
predict the output of forthcoming searches, the goal being to pro- 
vide information about the capabilities of future planet search 
instruments. With this mode, the flexibility of the code reaches 
its maximum, providing a wide range of possible applications. 

Once the synthetic planet population has been created, and 
assuming the characteristics of a given instrument, MESS_PM 
allows predicting the number of detections expected from a fu- 
ture facility. This provides informations on : 

1 . the expected frequency of planets 

2. the properties of these objects 

3. the kind of constraints that their observation can put on the 
planet formation theories. 

Furthermore, it also allows to test difi'erent instrumental config- 
urations and observational strategies that can be adopted, thus 
providing a tool to tune the instrument characteristics, in order 
to fulfil the requirements needed to access a certain domain in 
the parameter space, and reach the proposed science goals. 
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Table 2. Instruments for direct imagine of ex oplan ets considered in ou r 
analysis. References: BIO Beuzit et al/fcoiC); GOS'Green et al.VlOOS^; 
G07 Graham et al. (2007); KIO Kasoer et al. (2010); RIO Ri eke et al. 
dMoO ; SlO lStuiketalJ ( 120100 ; TlO lTrauger & Moodvi ( 120100 



Instrument 


Contrast 


Wavel. 


IWA 


Year 


Ref. 






(jim) 


(") 








S m ground-based 


telescopes 






VLT-SPHERE 


10-' 


0.9-1.7 


0.08 


2011 


BIO 


Gemini-GPI 










G07 


JWST 


NIRCAM 


10-^ 


2.1-4.6 


0.30 


2014 


G05 


MIRI 


10-* 


5-25 


0.35 




RIO 


1.5 m Space Coronagraphs 




10-" - 10-'" 


0.3-1.3 


0.08 


? 


TIO 


ELT's class instruments 


E-ELT-EPICS 


10-** - 10-" 


0.9-1.7 


0.03 


> 2018 


KIO 


E-ELT-METIS 


10-5 


2.5-20 


0.08 




SIO 



3.3.1. Comparison of future direct imaging instrument 
capabili 



ties 



As an example of the application of MESS_PM, we report the 
results of a comparison of the capability of a set of instruments 
for the direct imaging of exoplanets, planned for the next decade, 
which are briefly described in Tab. 2. 

Since the purpose of the presented analysis is purely illustra- 
tive, we adopted for each instrument an averaged detectability 
relation, taken from the reference indicated in Tab. 2, then us- 
ing only the ID approach. The sample of stars used is the one 
described in Sec. 12.2.11 and whose properties are summarized 
in Fig. [T] This sample was originally selected as a preliminary 
sample for the planet search survey to be done w ith SPHERE, 
the next generation planet finder of VLT dBeuzit et al.l2008l) . and 
it's therefore optimized for this kind of instruments, which possi- 
bly introduces some biases against some of the other instrument 
analyzed. The analysis was made using the reduced population, 
assuming five planet per star. 

The results of the analysis, showed in Fig. 12 and also sum- 
marized in Tab. 3, foresee that enormous progress that can be 
expected in the next decade. The available measurements are al- 
ready giving us indirect information on far away planets around 
young stars, but passing through the intermediate step of next 
generation image and finally with the advantage of ELT instru- 
ments we will have a wide view on planetary systems at different 
stages of their evolution. 



Table 3. Summary of expected detections from imagers in the next 
decade 



3.3.2. Foreseeing tlie synergies between (different 
tecliniques 

Once the RV and astrometric modules will be completed, 
MESS_PM will provide an estimate of both the direct and in- 
direct signatures of the presence of the planets, and thus be used 
to compare the outcomes of imaging with dynamical methods. 
These are interesting, because the latter allow determining the 
planet masses, thus eliminating the degeneracy with age, which 
is currently one of the major problems affecting direct detec- 
tions. Moreover, possible synergies between different discovery 
methods are becoming more and more likely, ELT's instruments 
representing the ideal link between direct and indirect detec- 
tions, covering both young, nearby systems discovered by next 
generation imagers and also meant to provide the first images of 
planets already detected by RV. 

Fig. 13 and 14 summarize the results of the preliminary ver- 
sion of the RV and astrometric modules of MESS_PM. The plan- 
ets showed are the same as in the lower right panel of Fig. 12. 

If confirmed, these results would suggest that the discovery 
space for EPICS at E-ELT overlaps well with those from ra- 
dial velocity (RV) instruments (HARPS at ESQ 3.6m telescope, 
ESPRESSO at VLT, and especially CODEX at E-ELT) as well 
as with that of GAIA dCasertano et al.ll2008l) . 

The RV module being still under test, and without having 
enough data to perform a consistent and accurate analysis of 
the performances and comparison between the instruments un- 
der scrutiny, this analysis is not meant to tell which instrument 
is going to provide the highest number of detection, but just at 
showing the potential of the further versions of the code. 
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Figure 13. Planets expected to be detected by EPICS (nearby 
sample) in the RV signal vs. period plane, compared with de- 
tection limits for RV instruments (HARPS, ESPRESSO and 
CODEX). The colour code is the same as in Fig. 1. 



Instrument Year Young Old Nept. Rocky 

Giants Giants 
Gr. based 8m 2011 tens few 
JWST 2014 tens few 

1.5m Space Coro. ? tens tens tens few 

ELT's > 2018 hundr. hundr. tens few 



3.4. Testing the influence of tine piiysicai inputs 

In this last section we present the results of some tests which goal 
was to show how MESS can be used to investigate the influence 
of the various physical parameters considered as inputs for the 
planet generation. In particular we focused on the eccentricity 
distribution and on the value of the planetary albedo. 
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Figure 12. Planets expected to be discovered by SPHERE (representative of planet finders on 8m class ground-based telescopes), 
JWST-MIRI, 1 .5m Space Coronagraphs, and EPICS/E-ELT (representative of 30-40m class telescopes) in the mass vs separation 
plane. Different colours are used for warm giant (orange), cold giant (red), Neptune-like (green), and rocky planets (blue) respec- 
tively. 



Eccentricity distribution 

Direct Imaging surveys are, by definition, mostly sensitive to 
planets in wide orbits. Also, planets on highly eccentric orbits 
could also be preferred targets, since they are more likely to 
be found farther out with respect to planets on a circular orbit 
with the same semi-major axis. This could led to a bias towards 
high eccentricity planets in our results. As mentioned in Sec. 2.2, 
the eccentricity distribution of the planets generated by MESS is 
uniform, and cut at e = 0.6. 

As a further check of the impact of the eccentricity on the 
DI results, we repeated the analysis of the TWA sample done in 
Sec. 3.2, by fixing the eccentricity for all the generated planets 
to a given value. The results are showed in Fig. 15. The black 
solid line shows the results of the standard setup, thus with the 
uniform eccentricity distribution cut at e = 0.6. The red, green, 
blue, purple and light blue lines show the outcomes of the sim- 
ulations done by fixing e = Q, e - 0.2, e - 0.4, e - 0.6 and 
e = 0.8, respectively. As expected, the higher eccentricity val- 
ues can lead to an higher fraction of detected planets, for a given 
semi-major axis value. 

This simple exercise shows not only that the standard setup 
of the MESS does not introduce any systematic bias towards 
high eccentricity, but also that the code allows us to easily take 



this kind of biases into account, if they are proven to be real, by 
changing the simulation parameters. 

As a final remark, it has to be said that the effect of the eccen- 
tricity is important only in the case of warm planets, as the ones 
that could be found around our TWA targets. As the age of the 
stars increases, the reflected light contribution to the planet con- 
trast becomes more and more important, thus counterbalancing 
the effect of the eccentricity. 

Albedo distribution 

As mentioned in Sec. 2.3, the Albedo of the planets in the syn- 
thetic population is randomly generated between 0.2 and 0.7. 
Especially for the cold planets, the value of the albedo can be 
a critical parameter for the planet detection. We then decided to 
perform a check to see how big is the impact on the simulation 
results. With an approach similar to the one used to test the ec- 
centricity effect (see Sec. 3.3.3), we performed different sets of 
simulation, being the only free parameter. We used an hypo- 
thetical G2V star (7 -, age - 4.5 Gyrs) at 20 parsecs as target, 
and the detection limits of EPICS (see Tab. 2). 

Fig. 16 shows the results of the standard setup {A^ randomly 
generated between 0.2 and 0.7, black solid line), together with 
the ones obtained by fixing A,i to 0.2, 0.35 (the Jupiter value. 
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Figure 14. Planets expected to be detected by EPICS (nearby 
sample) in the astrometric signal vs. period plane, compared 
with detection limits for astrometric satellites GAIA. Colour 
code is the same as in Fig. 1 . 
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Figure 15. Fraction of detected planet as a function of the semi- 
major axis value (AU) for different values of the planet eccen- 
tricity. The black solid line shows the results obtained with the 
standard setup. 



see Sec. 2.3.1), 0.5 and 0.7 (red, green, blue and purple line, 
respectively). 

As expected, the fraction of detectable planets is higher for 
higher albedo, all the other parameters being the same. 

This test confirms that the use of an uniform distribution al- 
low us not to favour planets with high albedo, the results being 
similar to the one obtained by using the Jupiter value. 

4. Summary and conclusions 

In this paper, we presented MESS (Multi-purpose Exoplanet 
Simulation System), a Monte Carlo tool for the statistical analy- 
sis and prediction of survey results for exoplanets. 

Our aim was to build an extremely versatile code, that could 
be used to test the outcomes of any instrument/technique for the 
detection of planets. We consider several assumptions on: 



Figure 16. Fraction of detected planet as a function of the semi- 
major axis value (AU) for different values of the planet albedo 
(Ai). The black solid line shows the results obtained with the 
standard setup. 

- The star population, and how to take into account the proper- 
ties of each star and their effect on either the characteristics 
of the planets or the instrument capabilities. The binarity as- 
pects is also included to take into account the possible effects 
of a stellar companion to the planet formation. 

- The planet population, providing the complete set of orbital 
elements and a large number of physical parameters of the 
planets (radius, temperature, luminosity, etc.), either gener- 
ated using the information coming from the analysis of the 
planets confirmed up to now (semi-empirical approach) or 
using the results of the planet formation theories (theoretical 
approach). 

- The predicted observables (luminosity and polarimetric con- 
trast, RV semi-amplitude, astrometric signal) 

- The synthesis of a planet population, that can be easily 
adapted to the purpose of the investigation 

- The final comparison with the detection limits, with the pos- 
sibility to combine the informations coming from different 
observing techniques, to select a sub-sample of detectable 
planets whose characteristics can then be investigated. 

The code is such that each and every one of these assumptions 
can be released and/or changed. This not only provides a tool 
which is independent from the models (e.g. the planet formation 
theory chosen if the theoretical approach is used, or the evolu- 
tionary models used to estimate the planet luminosity and radius) 
but also makes it relevant to test model prediction, as well as to 
constraint the properties of the known planets under different 
initial conditions. 

So far only the Direct Imaging module of the code has been 
extensively used, but the combination of various techniques is 
under test and will offer rich perspective for future combined 
studies of exoplanets. 

Three main applications of the MESS code have been shown: 

1. The Hybrid mode, built for the analysis of single objects, 
is presented in Sec l3.ll It can be used to probe the physical 
and orbital properties of a putative companion around a given 
system based on the combination of different techniques, and 
possibly a priori information on the orbit given the presence 
of other planets or of a circumstellar disk. 
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2. The SAM mode (Sec l3.2l l. optimized for the analysis of a 
large sample of stars, shows its full potential in Sec. 13.21 
by providing a detailed statistical analysis of a sample of 
stars observed with direct imaging. Both the agreement of 
the observations with the observed parameter distributions 
(Sec. l3.2.TT l and the planet formation theories (Sec. 13.2.3b are 
tested, using the semi-empirical and theoretical approach, re- 
spectively. 

3. The PM mode finally aims at the prediction of the outcomes 
of future searches, and can be used to tune not only the main 
instrument parameters, but even the observing strategy. 

However, an extensive use of the code requires a complete 
knowledge of the instrument under test, of all the error sources 
and of the detection capabilities. Then, to really extend the use 
of MESS to other facilities one should first properly set all the 
needed parameters. As already mentioned before, both the RV 
and the astrometric part are currently included in a very simplis- 
tic way. A better treatment of the dependence of the detectability 
with astrometry from the orbital parameters should be included. 
A rigorous treatment of the stellar jitter evaluation must be im- 
plemented to allow a better comparison between the imaging and 
radial velocity capabilities. Especially in the case of E-ELT in- 
struments, this would allow to better define the synergies be- 
tween the various channels, for a more focused observing strat- 
egy- 

Moreover, a precise measure of the stellar characteristics is 
also needed, in order to minimize the effects that errors on these 
parameters, such as the age of the system or the presence of stel- 
lar companions, can have on the analysis. 

Finally, the inclusion of an analysis of the planet stability in 
case of multiple objects is planned, together with an extensive 
use of the theoretica l approach, using the o utcomes of the most 
recent Bern models dMordasini et al.KOlOl) . 

Each technique performances vary with the star properties 
(age, mass, distance...), have different observables (luminosity, 
minimum mass, radius...), different observing strategies. It is 
therefore extremely important to take this into account to ring 
the maximum constrains, first on the properties of giant plan- 
ets (physical, orbital parameter space) that will actually entirely 
shape the planetary system architecture, then possibly on the 
telluric planets. A better characterization of the giant and tel- 
luric planet orbital and physical properties, including their de- 
pendency with the host properties, is critical for a better un- 
derstanding of their formation processes as various mechanisms 
may be at play (B oley 2009; Dodson-Robinson et al. 2009; 
IStamatellos & WhitworthI 12009). but also of their architecture 
and dynamical evolution. At the end, one additional and impor- 
tant issue is to understand the required physical conditions that 
will lead to the formation of telluric planets in habitable zone 
within planetary systems shaped by giant planets, and that wiU 
possibly lead to the formation of Life. 
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